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Рассматривается и решается задача получения аналитической зависимости выходной величины управляемо- 
го электрохимического сопротивления (УЭХС) от тока управления, а также выполнено математическое моде- 
лирование переходного процесса в выходной цепи УЭХС, содержащей линейную индуктивность или линей- 
ную электрическую ёмкость. 
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Введение. УЭХС представляет собой цилиндрический проводник длиною ( прямоугольного се- 
чения (подложка) ахг , изготовленный из стеклографита. Если на боковую поверхность УЭХС 
наносить электрохимическим способом тонкий слой металла, то сечение слоёного проводника 
будет увеличиваться, а его сопротивление уменьшаться. Ставится задача получить зависимость 
выходной величины УЭХС от тока (плотности тока) в цепи управления (рис. 1). 





Рис. 1. Схематическое устройство УЭХС 


Вывод зависимости А (Ё) от тока управления. Пусть р, и р, — удельные электрические 
сопротивления материалов подложки и наносимого осадка металла, $, и $,(х) — поперечные 
сечения подложки и осаждённого металла, АК, — элемент электрохимического сопротивления 


подложки длиною ах , аК, — элемент электрического сопротивления осаждённого металла дли- 
ною ах. Тогда 





р„ах 
ак, =——, 
- (1) 
р„ах 
ак, = ——-. 
Параллельное включение этих элементов приводит к такому эквиваленту 
В. . К 
а =". 
аК, + аК,„ (3) 
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Подставим (1) и (2) в (3) 


Р,@Х ХР, Ры 
аЕ = 5. $5„(х) ‘= 5 $„(х) ах = Р„Вы 
РобХ , Рыб р, Ры Ри5ы(Х) + Ры5% 


5, 5„() 5, 5„(х) 


Ш 








(4) 





Как видно из полученной формулы, дальнейшее исследование сводится к получению зависимости 
50) 

Пусть 0(х;Ё) — известное пространственно-временное распределение плотности элек- 
трического тока на управляющем электроде. Тогда 0(х;Ё).а.@х — ток управления в момент Е 


Е 

через площадку а-ах. | б(х;ЁК)а . ах .4Ё — заряд, прошедший через площадку а.аХ за отрезок 
0 

времени [0;# |. 


Если «К —  электрохимический эквивалент осаждаемого металла, то [1] 


Е 
ат =к [50а ах . Е — масса осаждённого металла на площадке а.ах , за отрезок времени 


[0;#]. | 


Е 
К! 5(х;Ба. ах. 
ат [54 ) 


Если у, — объёмная плотность осаждённого металла, то &М, = = — 
Ум Ум 


объём осаждённого металла на площадке а.4»Х за отрезок времени [0;2] . 





Считая форму осаждённого металла на поверхности бесконечно малой площади а.ах в 
виде параллепипеда, можно найти площадь поперечного сечения 5,(х) в точке х; 








Е 
ка| 5(х;Е)аЕ 
50 = ау, ‚ или $,(х) = — (5) 
ах У 
Подставим (5) в (4) 
ак = - РаРн ах. 
ка | 5(х; а (6) 
Ри +, 5 


Ум 
Выходная величина УЭХС легко может быть рассчитана с помощью интегрирования выражения 
(6) пох: 





6 

в(е) = | 1 ак. 
о ка[зоса (7) 
В 


Полученная формула (7) может быть использована при произвольном распределении плотности 
тока 5(х;Ё) в цепи управляющего электрода, в том числе и при равномерном, когда 


ы® 


о(Х;ЕЁ) = 
(х;6) те 


(8) 


где Г,(#) — ток управления. 
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Подставляя (8) в (7) получаем 





КЕ) ь РиРыё 


[% | Г(В (9) 


© 
Р» у, - 6 Ри?" 


м 


При постоянном токе управления Г, (Е) = Г, = сопЕ формула (9) ещё более упрощается: 





РоРий 
(Е = м . 
Е: (10) 
не 
Зависимость А (& в соответствии с формулой (10) принимает гиперболический характер (рис. 2). 
К 
ри 
5, 
0 Е 


Рис. 2. Зависимость К (& от времени при Г, (Е) = Г. = соп$ё 


Далее при расчётах переходных процессов в схемах, содержащих УЭХС, будем считать, 
что распределение плотности тока в цепи управляющего электрода является равномерным, т. е. 
будет использоваться формула (9). 


Математическое моделирование переходных процессов. Второй закон Кирхгофа для при- 
ведённой на рис. 3 схемы имеет вид 


аке) _ 
—_=Е» (11) 


Разрешая уравнение (11) относительно производной, получаем 
@_В-И(И(® (Е) +В) 


КЕВ(Е) + КОВ. +1 





р т : (12) 
По первому закону коммутации для данной цепи [2] 0-0) =/(0 +0) =0, т. е. 
(0) =0 (13) 
к(0 1 К 








И И ИИА 


























(ра 


ЦО 











Рис. 3. Включение А (#) с активно-индуктивной нагрузкой 


Задача коши (12), (13) легко может быть решена численно одним из известных методов, 
например, методом Эйлера. 


20 


Вестник ДГТУ. 2012. № 2 (63), вып. 2 








Введём в рассмотрение функцию 
Е -Г(Е) (В (Е) + А 
Г (Е (Е); (Е) = и (Е) ) 


Пусть Ё „ — время переходного процесса, п — количество дискретных временных точек 


(14) 





переходного процесса. Тогда временной шаг определяется соотношением 


В ие, (15) 
Г] 


Итерационная формула Эйлера для массива тока |1 } Я = 0, п] имеет вид 
а = + (8,51), (16) 
где К, =А (Ё,),/, =/(6,). 


Алгоритм расчёта переходного процесса в схеме, изображённой на рис. 3, можно предста- 
вить следующей блок-схемой (рис. 4). 


Ввод К, ; [, , параметров А (Ё), 


стандартной функции Г, (#), &„,М 






Вычисление /, (Е) 


Печать Ё;/ 


Вычисление интеграла 
Е 


Гру 





Вычисление А (Е) 





Останов 


Рис. 4. Алгоритм расчёта переходного процесса в цепи К (Е) -Ю-В 








При исходных данных Е =2В;р, =8.10?Ом.м;5, =10.10°м?;р„=1,72.103Ом.м; 
к =0, 6588 .10°кг / Кл; ( = 0,10м; 1, = 0,1А;у, = 8,9.103кг / мз; В, =5Ом;Ё =0,1гн была рассчи- 


тана зависимость / (Ё) ‚ представленная на рис. 5. 
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Рис. 5. Зависимость /(#) для схемы включения К (Ё)-К 4 Рис. 6. Включение А (#Ё) сактивно-ёмкостной нагрузкой 














Второй закон Кирхгофа для приведённой на рис. 6 схемы имеет вид 














Е 
КЕ + в, + [4 -Е,. (17) 
00 
Для получения дифференциального уравнения относительно / (Е) продифференцируем (17) по Ё. 
[6] .‚; (Е) 4 1. 
—_А(Ё Е К Е 
Е (В +К о ас ), 
откуда 
1 : К (Е) 
[6] С СЕ (18) 
=-КИ 
[4/8 К (Е) +К 
АК (ЕЁ) 
где ОЕ определяется при дифференцировании (9). 


Ввод К,; С, , параметров К (Ё), 


стандартной функции Г, (#), &„,М 











Вычисление /, (Ё) 


Печать Ё; / 


Вычисление интеграла 
Е 


Гр 


Вычисление 
к (2), (Е) 











С 


Рис. 7. Алгоритм расчёта переходного процесса в цепи к(-к -С 
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2 
р”рыК 
ое ой #(ё) 








Я >. (19) 
[м р, Е) +, >. 
У„/ 0 
По второму закону для схемы на рис. 6 можно записать 
и. (0-0)=И. (0+0) =0, 
откуда получаем при ё =(0+0) 
Ко) (0) +/ (0), =Е., 
т. е. 
. И Е 
16) к, о 


Таким образом, расчёт переходного процесса для / (Е) в схеме на рис. 6 сводится к решению за- 
дачи Коши (18), (20). 
Введём в рассмотрение функцию 
1 @8(Ё) 
© ЧЕ. (21) 


[веке Ея 


Итерационная формула Эйлера для массива тока |1, 1 =0, 0,7] имеет вид 








ра=б + (ВВГ, 
где 
ак (Е) 
СЕ 
Алгоритм расчёта переходного процесса в схеме, изображённой на рис. 6, можно представить 
следующей блок-схемой (рис. 7). 
При Е. =2В;К. =5Ом;С, =0,005ф была рассчитана зависимость тока /(Ё) (рис. 8) в схе- 





ве (Е = уве 


ме, представленной на рис. 6. 





Рис. 8. Зависимость тока /(Ё) для схемы включения К (Ё)- К, — С 
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Заключение. Полученная математическая модель УЭХС успешно может быть использована для 
расчёта переходных процессов в линейных и нелинейных схемах, содержащих этот элемент. В 
большинстве случаев, даже для линейных схем, эта задача решается только численно. 
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